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ABSTRACT 

We use gravitational lens models and X-ray spectral analysis of ten X-ray luminous galaxy clusters at z ~ 0.2 
to study the impact of cluster substructure on attempts to normalize the matter power spectrum. We estimate 
that unrelaxed clusters are 30% hotter than relaxed clusters causing tjg to be over-estimated by 20% if the cluster 
selection function is not accounted for correctly. This helps to explain the wide range in erg derived from differ¬ 
ent techniques, erg^ 0.6-1, and offers a physically motivated explanation for some of the discrepancy. We also 
identify two further systematics in our analysis: (i) extrapolation of small field-of-view mass measurements to 
the cluster virial radius and (ii) projection of 3-dimensional masses contained in numerical simulations to the 
2-dimensional information that is available from observations. We combine quantitative estimates of these two 
effects with our model fitting to estimate from the current data that as= 0.75 ± 0.05(statistical) ± 0.15(systematic), 
where the systematic error reflects the extrapolation and projection uncertainties. All three systematics (substruc¬ 
ture, extrapolation and projection) are fundamental to future cluster-based measurements of erg regardless of the 
techniques employed. However, we identify gravitational lensing as the tool of choice for such studies, because a 
combination of strong- and weak-lensing offers the most direct route to control the systematics and thus achieve 
an unbiased comparison between observation and theory. 

Subject headings: cosmology:observations — gravitational lensing — X-rays:galaxies:clusters — large-scale 
structure of the universe 

construct a robust cluster mass function from direct measure¬ 
ment. The local cluster X-ray temperature function, n(> T), 
has proved more accessible (e.g. Edge et al. 1990; Henry & Ar- 
naud 1991; Markevitch 1998; Blanchard et al. 2000; Pierpaoli 
et al. 2000; Ikebe et al. 2002). The X-ray temperature func¬ 
tion in conjunction with a robust mass-temperature calibration 
therefore offers an opportunity to constrain ug. 

Observational attempts to calibrate the mass-temperature re¬ 
lation typically rely on X-ray observations of clusters (e.g. 
Markevitch 1998; Nevalainen, Markevitch & Forman 2000; 
Finoguenov, Reiprich & Bohringer 2001; Allen, Schmidt & 
Fabian 2001, hereafter ASF; Reiprich & Bohringer 2002). De¬ 
spite the progress made by Allen (1998) in understanding X- 
ray based cluster mass measurements. X-ray techniques con¬ 
tinue to assume that all clusters are symmetric, equilibrium sys¬ 
tems. This is a major concern, because ~ 40-70% of galaxy 
clusters appear to be dynamically immature (e.g. Mohr et al. 
1995; Buote & Tsai 1996; Ota & Mitsuda 2002; Smith et al. 
2003, in prep., hereafter S03), and this immaturity has a mea¬ 
surable systematic impact on the normalization of the cluster 
mass-temperature and mass-luminosity relations (Ota & Mit¬ 
suda 2002; S03; see also Randall et al. 2002). 

In contrast, mass estimates based on gravitational lensing are 
insensitive to the physical nature and state of the cluster mass. 
Cluster lensing studies are therefore free from the symmetry 
and equilibrium assumptions that plague the X-ray studies. At¬ 
tempts to use lensing to calibrate the cluster mass-temperature 
relation have so far relied on previously published and/or crude 
cluster mass estimates (Hjorth et al. 1998; ASF; Viana, Nichol 
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1. INTRODUCTION 

The spectrum of cosmic matter fluctuations is an important 
constraint on theoretical models of structure formation (e.g. 
Press & Schechter 1974; Bond et al. 1991; Bower 1991; Kauff- 
mann & White 1993; Lacey & Cole 1993). The amplitude 
of the power spectrum is parametrized as erg, the linear-theory 
value of the rms fractional fluctuations in density averaged in 
spheres of 8 /j“*Mpc radius at z = 0. Several methods have 
been used to estimate ag: measurement of the abundance of 
galaxy clusters (Eke, Cole & Frenk 1996; Borgani et al. 2001; 
Reiprich & Bohringer 2001; Seljak 2001; Allen et al. 2002; 
Viana, Nichol & Liddle 2002; Pierpaoli et al. 2001), cosmic 
shear analyses (e.g. Bacon et al. 2002; Hoekstra et al. 2002; Re- 
fregier et al. 2002; van Waerbeke et al. 2002; Brown et al. 2002; 
Jarvis et al. 2002), cosmic microwave background (CMB) stud¬ 
ies (Sievers et al. 2002; Bond et al. 2002), combined analysis of 
galaxy redshift survey and CMB data (Lahav et al. 2002). Cur¬ 
rent estimates of ug range from ~ 0.6 to ~ 1.0, with claimed 
statistical uncertainties in the range Aug ~ 0.02-0.15. Overall, 
the situation is characterized by a lack of agreement between 
the results from different methods, or the same methods used 
on different samples, suggesting that systematic uncertainties 
probably lie at the heart of the current disagreement over the 
value of CTg. 

In this letter we investigate systematic biases in the use of 
cluster abundances to measure erg. In principal the mass func¬ 
tion of galaxy clusters, n(>M), should yield a direct constraint 
on CTg. However, it is not currently possible to measure clus¬ 
ter masses with the precision and in the numbers required to 
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& Liddle 2002). A major improvement on these pioneering 
studies would come from a precise and uniform analysis of a 
large objectively selected cluster sample for which high reso¬ 
lution space-based optical and X-ray data were available. In 
anticipation of such a program we conduct a pilot study us¬ 
ing S03’s Hubble Space Telescope (HST)/Chandra gravitational 
lensing survey of ten X-ray luminous galaxy clusters at z — 0.2. 
S03 made precise cluster mass and temperature measurements 
and thus constrained the high-mass end of the cluster mass- 
temperature relation. They also studied the dependence of this 
normalization on cluster substructure, concluding that unre¬ 
laxed clusters are, on average, 30% hotter than relaxed clusters. 
S03’s results therefore offer a unique opportunity to study the 
impact of cluster substructure on estimates of erg. 

We summarize S03’s results in §2, describe our model¬ 
ing and results in §3 and summarize our conclusions in §4. 
We express the Hubble parameter in terms of h, where Ho= 
100/zkms“^Mpc“*. We also adopt Ho= 0.3 and Ao= 0.7. 

2. HST/CHANDRA MASS-TEMPERATURE CALIBRATION 

S03 Studied a representative sample of ten of the most X-ray 
luminous clusters (Lx> 2 x lO'^h ^ ergs ', 0.1-2.4keV) in a 
narrow redshift slice at z = 0.21 ± 0.04, with line-of-sight red¬ 
dening of < 0.1 from the XBACs sample (Ebeling et 

al. 1996). Each cluster was typically observed for 3 orbits (i.e. 
7.5ks) through the E702W filter using the WEPC2 camera on¬ 
board HST. S03 used these data in conjunction with ground- 
based optical and near-infrared data (Smith et al. 2001, 2002), 
and the LENSTOOL software (Kneib 1993; Kneib et al. 1996; 
Smith 2002) to construct a detailed gravitational lens model of 
each cluster. 

Armed with these models, S03 measured M 2500 , the total pro¬ 
jected cluster mass within r 25 oo, i-C- the radius at which the 
density of matter in the clusters falls to p = P 2500 = 2500pc, 
where pc is the critical density required to close the universe.® 
S03 also used the models to divide their sample into relaxed 
(Afsub /Mtoi< 10%) and unrelaxed (Mjub /Aftot> 10%) clusters 
where Mtot is the total projected mass of the cluster within r 25 oo 
and Msub is the projected mass of the cluster within the same ra¬ 
dius that is not associated with the main centrally-located dark 
matter halo and cluster central galaxy. A complementary anal¬ 
ysis of archival Chandra and ASCA observations of eight and 
one of these clusters respectively also provided accurate mea¬ 
surements of the temperature of each cluster (Tx.tot) within a 
projected radius of r < 1 /i“'Mpc. We refer the reader to S03 
for further details of the modeling and analysis of these clusters. 

We plot S03’s mass and temperature measurements in Eig. 1. 
The open symbols show the individual clusters, and the filled 
symbols indicate the properties of the mean relaxed and unre¬ 
laxed cluster sub-samples. The mean temperatures of the re¬ 
laxed and unrelaxed clusters, are (Tx.tot) = 6.3 ±0.8 keV and 
(Tx.tot) = 9.2 ± 1.2 keV respectively, where the error bars are the 
uncertainties on the means and are the estimated from bootstrap 
re-sampling each sub-sample of clusters; we therefore expect 
these error bars to be over-estimates. The unrelaxed clusters 
appear to be systematically hotter than the relaxed clusters. 

Two of S03’s sample (A 383 - Smith et al. 2001; A 1835 
- e.g. Schmidt et al. 2001) have central cooling timescales of 
Aooi^ 10®years. This is in line with expectations from other 
representative samples of X-ray luminous clusters (Peres et al. 
1998). S03 therefore recalculated all of the cluster tempera¬ 


tures using an 0.05 < r < 1 /i“'Mpc annulus (i.e. excluding the 
cold core of the two extreme “cooling flow” systems). They 
found that, while the temperature difference is slightly reduced 
((Tx,ann)reiaxed = 6.9 ± 0.9 keV), it is robust to the exclusion of 
the central 50 /i“^kpc of each cluster from the temperature cal¬ 
culations. The 30% temperature difference therefore does re¬ 
flect a bona fide difference between the ambient temperatures 
of relaxed and unrelaxed clusters. We note that this affect is 
similar to a substructure-related bias found by Ota & Mitsuda 
(2002) in the cluster mass-luminosity relation. 


O 
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Fig. 1. — M2500, the projected mass within r25oo, versus the temper¬ 
ature of the intra-cluster medium for SOS’s sample of X-ray luminous 
clusters. The unrelaxed clusters are on average 30% hotter than the 
relaxed clusters, causing a previously unquantified structural bias in 
the normalization of the cluster mass-temperature relation. We also 
plot SOS’s and ASF’s mass-temperature relations, assuming a canoni¬ 
cal logarithmic slope of a = 2/3 for the SOS relation. The ASF relation 
agrees with the two cooling flow clusters in SOS’s sample (A 383 and 
A 1835 are the two open circles that lie within l-cr of the line). 

3. MODELING AND RESULTS 
3.1. Approach 

We construct a simple model to investigate the impact of 
S03’s results on estimates of erg. We begin by parameterizing 
the cluster mass-temperature relation: 

Txm (keV) = A(M25 oo/10'4/i-' Mq)“ (1) 

where Tx^ot and M 2500 are as defined in §2, and A,a are the 
normalization and logarithmic slope respectively. We first con¬ 
vert the Jenkins et al. (2001) mass function to a temperature 
function. This conversion includes the following elements: a 
mass-dependent concentration index (Eke, Navarro & Stein- 
metz 2001), conversion of three-dimensional masses from the 
simulations to projected two-dimensional masses (Hjorth et 
al. 1998) and an observational mass-temperature normalization 
(S03, ASE). We then fit this model temperature function to the 
observed temperature function (Edge et al. 1990) using a single 
free parameter, erg. 

Our model also contains the following parameters: {Oq, Aq, 
r, aj } where Oq and Aq are the matter and vacuum energy 


'At z = 0.2, r 25 oo corresponds to the edge of the /fiS7’AVFPC2 field of view for the most massive clusters in S03’s sample. 
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densities of the Universe at z = 0, F is the spectral shape param¬ 
eter for the power spectrum and aj is the scatter in log(7’x^tot)- 
We focus our attention on the dependence of a& on A and to a 
lesser extent on ctj; we therefore adopt “concordance” values 
for the remaining parameters: Uo= 0.3, Ao= 0.7, F = 0.2 (e.g. 
Efstathiou et al. 2002), a = 2/3 (e.g. ASF). We stress that we 
adopt a fixed value of r2o= 0.3, and therefore do not investigate 
the CTg-Uo degeneracy. 

3.2. Model Fitting 

We use two independent mass-temperature calibrations to 
normalize our models. We begin with S03’s normalization, and 
adopt the values of A and ar relevant to their entire sample: 
A = 4.4, CTr = 0.1 (see solid line in Fig. 1). We compute a model 
temperature function and fit it to the observed temperature func¬ 
tion (Fdge et al. 1990), obtaining a best-fit of ag,= 0.75 ±0.05 
where the uncertainty is the statistical uncertainty on the fit. We 
plot this best-fit model and the observed temperature function in 
Fig. 2. Next, we turn to ASF’s cooling flow mass-temperature 
relation. These authors observed a sample of seven cooling 
flow clusters with Chandra, and used these data to normalize 
the mass-temperature relation. We convert ASF’s cooling flow 
mass-temperature relation into the form required for our model: 
A = 2.6, ar — 0.03. Using these values, we construct a model 
temperature function and fit it to the observed function, obtain¬ 
ing (78= 0.91 ±0.07. This model (Fig. 2) fits the data less well 
than the S03-based model, with the largest residuals occurring 
at high temperatures. 



kT(keV) 

Fig. 2. — We plot the Edge et al. (1990) cluster temperature function 
for both all clusters and cooling flow clusters (defined as containing 
a line-emitting central galaxy), together with the best-fit model tem¬ 
perature functions that are normalized with the S03 and ASF mass- 
temperature relations. When a cooling flow cluster based normal¬ 
ization is applied to a representative sample of clusters, ag is over¬ 
estimated by ~ 20%. However, when the cluster selection function is 
accounted for properly in both the model normalization and the ob¬ 
served temperature function, consistent values of ag are obtained. 

A simple interpretation of these two models is that it is im¬ 
portant to understand the cluster selection function when using 
cluster abundances to measure a&. Specifically, using a cooling- 
flow cluster mass-temperature normalization when studying the 


temperature function of a representative sample of clusters may 
cause (Tg to be over-estimated by approximately 20%. 

We test this interpretation by fitting the cooling flow normal¬ 
ized model to a temperature function that describes just cooling 
flow clusters. We first use the observed correlation between 
line emission from cluster central galaxies and short cooling 
timescales (fcooi^ 10^ years - Fdge et al. 1992; Peres et al. 
1998; Crawford et al. 1999) to construct a “cooling flow only” 
temperature function from the Fdge et al. (1990) sample. We 
then fit the cooling flow model to the cooling flow data and 
obtain a best-fit value of as= 0.74 ± 0.05, which agrees with 
the value obtained from the original model that was normalized 
with S03’s results. We plot this best-fit model and the relevant 
data in Fig. 2. This model confirms our interpretation that clus¬ 
ter substructure is an important and previously unidentified sys¬ 
tematic effect at the 20% level when using cluster abundances 
to constrain erg. 

3.3. Further Uncertainties 

The low value of erg obtained in §3.2 is similar to a number 
of other recent results that favor (Tg~ 0.6-0.8 (e.g. Seljak 2001; 
Reiprich & Bohringer 2001; Borgani et al. 2001; Allen et al. 
2002; Viana, Nichol & Fiddle 2002; Lahav et al. 2002; Brown 
et al. 2002; Jarvis et al. 2002). However, several other uncer¬ 
tainties need to be investigated before a reliable conclusion on 
the value of erg from cluster abundance determinations can be 
drawn. 

Firstly, we highlight the extrapolation of S03’s lens models 
from 72500 (i e. approximately the edge of the //57’/WFPC2 field 
of view at z ~ 0.2) to the cluster virial radii as a key systematic 
uncertainty in our analysis. Bardeau et al. (in prep.) investi¬ 
gate this effect in detail through their weak-shear analysis of 
panoramic (28' x 42') CFH12k BRI-hand imaging of the S03 
cluster sample. Prior to the completion of this wide-held anal¬ 
ysis, we note that weak lensing analyses of individual clusters 
(e.g. King, Clowe & Schneider 2002) are unable to discrimi¬ 
nate between isothermal (p (x r~^) and Navarro, Frenk & White 
(1997) (p (X r~^) profiles on large scales. We therefore exploit 
this lack of discriminatory power to make a conservative es¬ 
timate of this systematic uncertainty. We integrate both pro¬ 
files over the radial range 0.25 ^ r ^ 1.5 fi“^Mpc (i.e. the dy¬ 
namic range over which we are extrapolating), and estimate that 
the uncertainty in profile shape introduces an uncertainty in the 
virial mass estimate for an individual cluster of ~ 30%, which 
translates into an uncertainty in cluster temperature (assuming 
of ~ 20%. This equates to a ~ 10% “extrapolation” 
systematic uncertainty in erg. 

Secondly, we identify the projection of three-dimensional 
cluster masses from numerical simulations to observed two- 
dimensional masses (§3.1) as a further source of systematic un¬ 
certainty. As Hjorth et al. (1998) discuss, the magnitude of this 
uncertainty depends on the slope of the cluster density profile 
at small radii. Recent observational results (Smith et al. 2001; 
Sand, Treu & Fills 2002; Dahle, Hannestad & Sommer-Farsen 
2002) indicate that there may be substantial intrinsic scatter in 
this slope, appearing to contradict theoretical claims for a uni¬ 
versal profile (e.g. Navarro, Frenk & White 1997). Given these 
complications, we conservatively adopt a further 10% “projec¬ 
tion” systematic uncertainty in erg. 

In summary, although S03’s detailed lens models allow the 
“substructure” systematic to be accounted for properly and (to 
first order) eliminated from our analysis, “extrapolation” and 
“projection” uncertainties combine to produce a ~ 20% sys- 
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tematic that we are unable to control with the current dataset. 

4. SUMMARY AND DISCUSSION 

We have used S03’s substructure-dependent cluster mass- 
temperature normalization to investigate the impact of cluster 
substructure on estimates of cts. We find that when a cooling 
flow cluster mass-temperature normalization is applied to the 
general cluster population, a& is over-estimated by 20%. A 
clear understanding of the cluster selection function is there¬ 
fore fundamental to attempts to constrain ag with cluster abun¬ 
dances. The simple X-ray luminosity-limited selection of S03’s 
sample (§2) enable us to account for this “substructure” sys¬ 
tematic from our analysis and thus to estimate that ag= 0.75 ± 
0.05(statistical). However, before we conclude that ag= 0.75, 
we highlight two further systematic effects which may bias our 
analysis: extrapolation of S03’s small field-of-view lens models 
out to the cluster virial radii, and uncertainties in the relation¬ 
ship between three-dimensional mass information contained in 
numerical simulations and the two-dimensional mass informa¬ 
tion that is available from observations. We estimate conser¬ 
vatively that these effects combine to produce a further 20% 
systematic uncertainty, and therefore we conclude from the 
present datathat erg = 0.75±0.05 (statistical) ±0.15 (systematic). 
We also note that the recently reported discrepancies between 
XMM- and Chandra-hased cluster temperature measurements 
(Schmidt et al. 2001; Majerowicz et al. 2002; Markevitch 2002) 
may introduce further uncertainty into cluster abundance deter¬ 
minations of ag. 

Our 20% “substructure” systematic is similar to the discrep¬ 
ancy between the canonical value of ag'^ 0.9-1 (e.g. Eke, Cole 
& Frenk 1996; Pierpaoli et al. 2001; Bacon et al. 2002; Bond 
et al. 2002; Hoekstra et al. 2002; Refregier et al. 2002; van 
Waerbeke et al. 2002) and recent claims for ug'^ 0.6-0.8 (Sel- 
jak 2001; Reiprich & Bohringer 2001; Borgani et al. 2001; 
Allen et al. 2002; Brown et al. 2002; Jarvis et al. 2002; La- 


hav et al. 2002; Schuecker et al. 2002; Viana, Nichol & Lid- 
die 2002). Our results therefore offer a physically motivated 
explanation for some of this discrepancy. Independent confir¬ 
mation of this comes from Randall et al.’s (2002) semi-analytic 
study of the effect of cluster mergers on the observed luminos¬ 
ity and temperature functions, and thus on the inferred clus¬ 
ter mass function. Randall et al. predict that cluster mergers 
boost the observed temperature function and can cause ag to be 
over-estimated by 20% if hydrostatic equilibrium is assumed 
for non-equilibrium clusters, in agreement with our observa¬ 
tional results. 

All three systematics discussed in this letter (substruc¬ 
ture, extrapolation and projection) affect the ability of cluster 
abundance techniques to measure ag accurately, regardless of 
whether gravitational lensing or X-ray techniques are used to 
measure the cluster masses. However, the insensitivity of grav¬ 
itational lensing to the physical nature and state of the cluster 
matter means that a combined strong- and weak-lensing study 
of a large, objectively selected sample of clusters should be the 
tool of choice for future cluster abundance studies. 
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